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ABSTRACT 

We present Chandra observations of 12 galaxies that contain supermassive black holes with dynamical mass 
measurements. Each galaxy was observed for 30 ksec and resulted in a total of 68 point source detections in the 
target galaxies including supermassive black hole sources, ultraluminous X-ray sources, and extragalactic X- 
ray binaries. Based on our fits of the X-ray spectra, we report fluxes, luminosities, Eddington ratios, and slope 
of the power-law spectrum. Normalized to the Eddington luminosity, the 2-10 keV band X-ray luminosities of 
the SMBH sources range from 10"^ to 10"^ and the power-law slopes are centered at ^ 2 with a slight trend 
towards steeper (softer) slopes at smaller Eddington fractions, implying a change in the physical processes 
responsible for their emission at low accretion rates. We find 20 ULX candidates, of which six are likely 
(> 90% chance) to be true ULXs. The most promising ULX candidate has an isotropic luminosity in the 
0.3-10 keV band of LO+Jj;^ x lO"*" erg s"'. 

Subject headings: accretion — accretion disks — black hole physics — galaxies: general — galaxies:nuclei 
— X-rays: galaxies — X-rays: general 



1. INTRODUCTION 

The prevalence of large black holes at the centers of galax- 
ies (e.g., Richstone et al. 1998) and their role as the central 
engines of active galactic nuclei (AGNs; e.g., Rees 1984) 
has been well appreciated. The possibility of convolution 
of black holes and their host galaxies, particularly through 
self-regulated growth and feedback from accretion-powered 
outflows (Silk & Rees 1998; Fabian 1999), has focused re- 
cent observational, theoretical, and computational effort (e.g., 
Schawinski et al. 2007; Di Matteo et al. 2005). 

The detailed microphysics at play in accretion-powered 
feedback, however, is not yet understood. Accretion onto a 
black hole is thought to proceed via an accretion disk, which 
launches relativistic jets and outflows from its inner regions 
(e.g., Lynden-Bell 1978; Blandford & Payne 1982; Blandford 
& Znajek 1977). In the same way that accretion disk prop- 
erties appear mostly to scale with black hole mass, so do the 
length- and time-scales of jets from stellar-mass black holes 
to supermassive black holes (SMBHs). 

The radio emission, X-ray emission, and mass of an accret- 
ing black hole are empirically related through what is some- 
times called "the fundamental plane of black hole accretion" 
(Merloni et al. 2003; Falcke et al. 2004). Radio emission 
coming from synchrotron emission in the jets clearly must 
depend on the amount of matter accreted towards the black 
hole, at which point it is turned into an outflow. X-ray emis- 
sion has several possible origins, including from the accretion 
disk corona, or from jets. The X-ray emission depends on the 
accretion rate and the source compactness, itself a function of 
the size of and thus mass of the black hole. Thus, it is natural 
to expect some form of mutual covariance among these three 
quantities. The relatively small scatter in the relation, which 
spans over 8 orders of magnitude in black hole mass, how- 
ever, suggests that accretion and outflows are self-regulated 



in a similar way across in all black holes. 

Since the discovery of the fundamental plane (Merloni et al. 
2003; Falcke et al. 2004), there has been a concerted effort to 
understand it, both from an observational perspective, which 
focuses on the universality and extent of the relation, and 
also from a theoretical perspective, which has focused on un- 
derstanding the mechanisms of jet production, constraining 
which radiative mechanism (or mechanisms) drive the corre- 
lation and at what efficiency for jets and accretion inflow pro- 
cesses (Merloni et al. 2006; Kording et al. 2006; Wang et al. 
2006; Li et al. 2008; Yuan et al. 2009; Plotkin et al. 201 1; de 
Gasperin et al. 2011). 

In Gultekin et al. (2009a), we performed an archival Chan- 
dra analysis of all SMBHs with primary, direct black-hole- 
mass measurements and available radio data. We make a dis- 
tinction between mass measurements that are "primary and 
direct" such as stellar dynamical (e.g., Gultekin et al. 2009b), 
gas dynamical (e.g., Barth et al. 2001), and megamaser mea- 
surements (e.g., Kuo et al. 2011) and reverberation mapping 
measurements (e.g., Bentz et al. 2006). The reverberation 
mapping technique is direct in that it directly probes a black 
hole's surrounding environment but is secondary because it 
must be normalized to the quiescent population through host- 
galaxy scaling relations (Onken et al. 2004). There are sev- 
eral advantages to using only sources with primary black hole 
masses. First, statistical techniques that make use of mea- 
surement errors can be faithfully used with the actual errors 
on black hole mass rather than on an inferred scatter to a host- 
galaxy scaling relation (e.g., Gultekin et al. 2009c). Second, 
by focusing on black holes with known, dynamical masses, 
we can measure true Eddington fractions and understand how 
accretion processes depend on this. Third, and perhaps most 
importantly, the results can be used to directly calibrate an es- 
timator for black hole mass based on nuclear X-ray and radio 
measurements. 
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In this paper we make the first step in completing the X- 
ray and radio survey of black holes with primary, direct mass 
measurements with 12 new Chandra observations. This paper 
will be complimented with one using Extended Very Large 
Array (EVLA) observations of the same sources. While the 
sample was designed to study the fundamental plane, a Chan- 
dra survey of SMBHs with dynamical mass measurements 
provides an interesting look at the X-ray properties of low lu- 
minosity AGNs (LLAGNs). LLAGNs do not appear to be 
scaled-down versions of more luminous Seyferts and quasars, 
but instead display very different accretion physics (Ho 2008). 
X-ray emission is generally considered a probe of accretion 
power, wherever it is ultimately produced. Because of the 
expected differences in physics at low accretion rates. X-ray 
studies of LLAGNs should provide as direct a look as possi- 
ble. 

In this and other X-ray observational work on LLAGNs, a 
key measure of the difference in their energetics is the how 
spectral shape changes with accretion rate. One simple, di- 
rect, and purely observational diagnostic is the change in 
slope of hard X-ray emission (2-10 keV) with hard X-ray Ed- 
dington fraction. This has been examined in both X-ray bi- 
naries (Corbel et al. 2006) and AGNs (Shemmer et al. 2008; 
Gu & Cao 2009; Winter et al. 2009; Constantin et al. 2009; 
Younes et al. 2011). One potential limitation of the AGN 
studies is the use of secondary mass estimates in determin- 
ing Eddington fraction. With our sample, we are not limited 
by this, and we are able to test down to very low accretion 
rates. 

In Section 2 we describe the sample selection, observations, 
data reduction, and spectral fitting. In Section 3 we detail the 
results of the SMBH sources as well as serendipitous detec- 
tion of ultra luminous X-ray sources (ULXs) and other point 
sources in the target galaxies. 

2. OBSERVATIONS AND SPECTRAL FITTING 

In this section we describe our experimental method from 
sample selection, though data reduction, point source detec- 
tion, and spectral extraction and fitting. All data reduction 
was done with CIAO version 4.3 and calibration databases 
(CALDB) version 4.4.3, and spectral fitting was done with 
XSPEC version 12 (Arnaud 1996). Reduction was done with 
the distributed Level 2 event files, which were processed at 
different times but all between 2010 Apr 14 (NGC 4486A) 
and 2010 Dec 12 (NGC 4291). 

2.1. Sample Selection 

In Giiltekin et al. (2009a), we analyzed archival Chandra 
data of all SMBHs with dynamical masses. The parent sam- 
ple was the list of "secure" black hole mass detections in 
Giiltekin et al. (2009c). Of the approximately 50 black holes 
in that list, 21 had only limits on their X-ray luminosity or 
else were not observed with Chandra at all. Of these 21, we 
identified 13 whose nuclear flux could potentially be deter- 
mined with a 30 ks Chandra observation. The others had 
large amounts of contaminating hot gas or existing archival 
data indicated that they were too faint to be observed with 
30 ks exposures. We were granted 12 of these observations 
with joint Extended Very Large Array time. The rest should 
be observed with X-ray Multi-mirror Mission-Newton, which 
has a larger effective area, or are extremely faint objects that 
have not been detected with deep Chandra exposures. Ta- 
ble 1 summarizes target galaxies and masses of the central 



Table 1 

Chandra Observations 



Obs ID 




l-rn 1 Q 

yja.la.Ky 


Dist 


Ret" 


^BH 


Ref 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


11775 


30.04 


NGC 1300 


20.1 


Dl 


7.1^1^ X 10' 


Ml 


2.53 


11776 


30.05 


NGC 2748 


24.9 


Dl 


A.l+\l X 10' 


Ml 


1.55 


1 1777 


29.55 


NGC 9778 


24 2 


D2 


1 6+"-'' X 10' 

10 


M2 


1.61 


11782 


29.04 


NGC 3384 


11.7 


D2 


1.8+Ui X 10' 


M2 


2.94 


11778 


30.16 


NGC 4291 


25.0 


D2 


3.2+0:5 X 10** 


M2 


2.87 


11784 


30.18 


NGC 4459 


17.0 


D3 


7.4+1:^ X 10' 


M3 


2.67 


11783 


29.05 


NGC 4486A 


17.0 


D3 


X 10' 


M4 


2.00 


11785 


31.38 


NGC 4596 


18.0 


D4 


8.4+^-° X 10' 


M3 


1.43 


11779 


33.08 


NGC 4742 


16.4 


D2 


1.4+;j:^xl0' 


M5 


3.43 


11780 


29.05 


NGC 5077 


44.9 


D5 


X 10** 


M2 


3.05 


11781 


30.05 


NGC 5576 


27.1 


D2 


1.8^U;3^10** 


M7 


2.44 


11786 


29.04 


NGC 7457 


14.0 


D2 


4.1+[-^ X 10'' 


M2 


4.60 



References. — (Dl) Atkinson et al. 2005, (D2) Tonry et al. 2001, (D3) 
Mei et al. 2007 adopted to 17.0 Mpc, (D4) Tully 1988, (D5) Faber et al. 1989 
group distance, (Ml) Atkinson et al. 2005, (M2) Gebhardt et al. 2003, (M3) Sarzi 
et al. 2001, (M4) Nowak et al. 2007, (M5) Giiltekin et al. in preparation, (M6) de 
Francesco et al. 2008, (M7) Giiltekin et al. 2009b. 



Note. — We summarize the observations and the targets. The columns are (1) 
the Chandra observation identification number; (2) the exposure time in units of 
ks; (3) the targeted galaxy; (4) the distance of the galaxy in units of Mpc; (5) a 
reference code for the distance; (6) the mass of the central black hole in units of 
Mv, with IfJ uncertainties; (7) a reference code for the mass measurement; (8) the 
Galactic column towards each source in units of 10-" cm^^ (Bajaja et al. 2005). 
The distances have been scaled from the indicated literature values to a common 
Hubble parameter of h = 0.7, and NGC 4459 and NGC4486A have been set at a 
Virgo distance of 17.0 Mpc. 

black holes. All but one galaxy are within 30 Mpc, and the 
black hole masses span a wide range: 4.1 x 10^-8.0 x 10** M,,.,. 
The only galaxy with a nuclear activity classification is NGC 
5576, which is classified as a Low Ionization Narrow Emis- 
sion Region (LINER) with broad Balmer lines (Ho et al. 1997; 
Veron-Cetty & Veron 2006). 

One potential selection effect that arises from requiring a 
dynamical mass measurement for inclusion in our sample is 
that it is biased to low Eddington rates. In general, the meth- 
ods used to measure black holes in this sample are hampered 
by contamination from a strong AGN contribution in the same 
bands used for the mass measurement. For this reason, the re- 
sults and conclusions we draw from this study only apply to 
the very low Eddington rates (~ 10^^-10^^) that our data ad- 
equately cover 

2.2. Point source detection 

To detect points sources in the field of each galaxy, we used 
the wavdetect tool on the whole image in the whole Chan- 
dra band at full resolution, run with the large_detect.pl wrap- 
per script provided by the ACIS team^ The wrapper script 
splits the image into multiple, overlapping sub-images and 
runs wavdetect^ on each region, identifying multiply detected 
sources in the overlapping portions to minimize edge effects 
in the wavlet algorithm. Wavdetect was used default settings 
for detection and background rejection thresholds and with 
searches in wavelet radii (i.e., the "scales" parameter) of 1, 2, 
4, 8, and 16 pixels. The whole image was used rather than 
just the region near the galaxy because 

' See http;//goo.gl/hk7e6 

^ It also runs celldetect, but we only use the wavdetect result. 
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For each Chandra image we considered only the point 
sources that were located on the galaxy as determined by 
comparison to 2MASS and DSS images of the galaxy. 
The size of the galaxy was determined manually for each 
galaxy by looking at the images and is roughly compara- 
ble to using isophotes of 25 mag arcsec^^ in the B-band or 
20 mag arcsec^^ in the K band. The range of values for di- 
ameters at a particular surface brightness in a particular band 
for each galaxy could vary by as much as 10%, motivating 
our choice for a simple, manual estimate. In total, we found 
68 point sources over the 12 images. We list all point sources 
detected in table 5 with a running identification number, lAU 
approved source name, J2000 coordinates, and net count rates 
in the 0.3-1, 1-2, and 2-10 keV bands. We show positions of 
all sources in the Chandra images as well as locations on the 
DSS images of the galaxy in Figure 1. 

2.3. Astrometry 

Because we are principally interested in the nuclear X-ray 
source of each galaxy, we have taken extra care in identifying 
which point source is the nuclear source. For most galaxies 
it was unambiguous: there was only one point source consis- 
tent with the optical or infrared center of the galaxy. In three 
galaxies (NGC 4486A, NGC 4596, and NGC 4742) there 
were two X-ray point sources that were consistent with the po- 
sition of the galaxy's center. In a fourth galaxy (NGC 3384), 
there was faint emission that was not detected by wavdetect 
as a point source, but which we considered for potential con- 
fusion. Typical separation of the confusing sources was 1 to 
3", which requires only modest absolute astrometric correc- 
tions. For the four galaxies, we registered the Chandra im- 
ages to Sloan Digital Sky Survey (SDSS) or Deep Near In- 
frared Survey (DENIS) coordinates using background AGNs 
that appeared in both the X-ray and optical/infrared images. 
Each image had between 3 and 6 sources for registration. The 
optical/infrared center was obvious in each image so that we 
could unambiguously determine which source was closer to 
the galaxy center In the end, we were able to determine that 
the Chandra coordinates were sufficient for determining the 
nuclear sources. The nuclear source in each galaxy is identi- 
fied in Table 5 in the classification column as "Nuc." 

One possible source of confusion is if an X-ray binary was 
located at the position of the nuclear center. To estimate this 
amount of contamination, we adopt the Grimm et al. (2003) 
universal cumulative luminosity function of high-mass X-ray 
binaries (HMXBs) in a given galaxy as a function of star for- 
mation rate, SFR: 



A^(>L) =5.4(SFR/Me yr"^) (L/10^^ erg s"^' 



-0.61 



0.038 
(1) 

where the second term in the brackets comes from assuming 
an upper limit of Lx =2.1 x 10^'' erg s"^ for a HMXB and 
can be neglected for our calculations. We take the lumin- 
soity of each nuclear source; e.g., for the dimmest source, 
Lx = 3 X 10^^ erg s"' so that N{L > 3 x 10^^ erg s"') = 
2.8(SFR/M,:, yr^'). We calculated the star formation rate, 
SFR, assuming the (Rosa-Gonzalez et al. 2002) estimator. 



SFR = 4.5 X lO"^'* Ma yr"^ (LpiR/erg s" 



(2) 



where the far-infrared luminosity is calculated as LpiR = 
1 -iXFx at A = 60 /im. This estimator is based on the Ken- 
nicutt (1998) prescritpion and agrees with other, similar pre- 



scriptions (e.g.. Bell 2003). Using far-infrared may overesti- 
mate the star-formation rate in early-type galaxies but is con- 
servative for our argument. We based our calculations on 
the the IRAS 60 jUm flux density (Knapp et al. 1989) for aU 
galaxies except NGC 5077 for which we substituted the MIPS 
70 /im flux density (Temi et al. 2009). The galaxy with the 
highest star-formation rate in our sample is NGC 2748 with 
SFR = 2.0 yr- ' , and all but 4 have SFR < 0. 1 M, yr" ' . 
If we assume that the X-ray binaries follow the light in the 
galaxy, then we can estimate the number of potential contam- 
inating sources. This is an overestimate of central contamina- 
tion probability for galaxies like NGC 1300 for which X-ray 
binaries follow the high rate of star formation in the spiral 
arms. The starlight in the central 1", roughly our astrometric 
uncertainty, ranges from approximately 10% for NGC 1300 to 
less than 1%. Assuming 10% for all galaxies, we can calcu- 
late the expected number of contaminating sources at the cen- 
ter each galaxy. The values range from 0.04 for NGC 2748 
to less than 0.01 for the majority. These numbers are con- 
servative because the actual amount of light at the center of 
each galaxy is less 10%, and the star formation rate may be 
over-estimated in the early-type galaxies. 

For early-type galaxies, a greater concern is the chance po- 
sitioning of a low-mass X-ray binary (LMXB), which gener- 
ally traces an older population (Kim & Fabbiano 2004). In 
this case, we use the results of Kim & Fabbiano (2004), who 
found that the cumulative number of LMXBs greater than a 
luminosity, Lx to scale as N{> Lx) ^ with a total lumi- 
nosity of the LMXB population scaling with the galaxy lumi- 
nosity Lk as Lpop = (0.2±0.08) x 10^° erg s-^Lk/Lk,. This 
is the most concern for NGC 5077, which has the largest total 
K-h&nA luminosity of L^ = 8.5 x lO'" L/f v^, (Skrutskie et al. 
2006), and for NGC 4486A, which has the lowest nuclear lu- 
minosity (Lx = 2.5 X 10^^ erg s^^) and a /T-band luminosity 
of La: = 1.5 X 10^" Lk., (Skrutskie et al. 2006). Again as- 
suming that LMXBs follow the light, the expected number of 
LMXBs that would be as bright or brighter than the identified 
nuclear sources is less than 0.005 and 0.04 for NGC 5077 and 
NGC 4486 A, respectively. It is lower for all other galaxies. 
Because each galaxy type is predominantly affected by con- 
tamination from only one of HMXBs or LMXBs and not both, 
the combined expected contamination by both are then all less 
than 0.05. Thus we expect not to have incorrectly identified 
an X-ray binary as an SMBH source. 

2.4. Spectral reduction 

We followed the standard pipeline in reduction of all data 
sets, using the most recent Chandra data reduction soft- 
ware package (CIAO version 4.3) and calibration databases 
(CALDB version 4.4.3). There was no significant background 
flaring, so there was no need for filtering. We used the CIAO 
tool psextract to extract the point-source spectra. Since our 
observations all used the Advanced CCD Imaging Spectrom- 
eter (ACIS), we ran psextract with the mkacisrmf tool to cre- 
ate the response matrix file (RMF) and with mkarf set for 
ACIS ancillary response file (ARE) creation. Source regions 
were circles with centers at the coordinates given by wavde- 
tect. The radii of the regions were the semimajor axes of the 
wavdetect error ellipses. This accounts for two effects: (1) 
the uncertainty in the location of the source and (2) the degra- 
dation in the point spread function off axis. For background 
regions, we used annuli with inner radii just larger than the 
source region radius and outer radii that was typically 14 pix- 
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Figure 1. Each pair of plots sliows tlie Chandra image (left) with detected point sources labeled with their running identification number and the DSS image 
(right) with locations of Chandra point sources labeled. The circles are centered at the location of the extraction region and are drawn with an arbitrarily sized 
radius for clarity rather than indicating the size of the extraction region. The bar under the galaxy name is 30" long. 
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els larger 

2.5. Spectral fitting 

We modeled the reduced spectra using XSPEC12 (Amaud 
1996). For some sources we used statistics and for others 
C-stat statistics (Cash 1979). The decision on which statistics 
to use was based on the number of photons available for bin- 
ning in energy bands. If binning the spectra in energy so that 
each bin contained a minimum of 20 counts resulted in five 
or more bins, we used statistics; otherwise we used C-stat 
statistics with unbinned spectra. There is necessarily a loss of 
information when binning the spectra, but we found when us- 
ing both types of statistics on sources that had enough counts 
to support it, the resulting parameter estimates were always 
consistent within la. We report the results from x^ statistics 
on account of the intuitive nature of the goodness-of-fit with 
X^ statistics. 

All spectra were modeled with a photoabsorbed power-law 
model with the intrinsic flux modeled directly as one of the 
parameters, i.e., the XSPEC model used was phabs(cflux 
* powerlaw), with the cf lux component normalized to the 
2-10 keV band. The spectra were fitted from 0.3 to 10 keV. 
While most of the photon counts above ^ 7 keV are proba- 
bly dominated by background, the fitting methodology takes 
this into account. We tested our results for sensitivity to the 
adopted upper energy cutoff and found that as long it was 
greater than 5 keV, our results were robust in the sense that 
they changed by far less than the Id uncertainties. When fit- 
ting the spectrum, we set a hard minimum for the absorption 
to be the Galactic value for A^h towards each source (Bajaja 
et al. 2005), but in calculating the uncertainty in the column, 
we allowed it to drop below this value. The results of the spec- 
tral fits are presented in Table 2 with best-fit parameters and 
la (68%) uncertainties for A^h, logf2_io, and F (the power- 
law index). Some sources did not have sufficient counts to 
produce reliable fits and are not included in the table. Fits to 
three sources (5, 14, and 62) were unconstrained on at least 
the spectral index parameter, F, and we consider these fits ap- 
proximate. 

3. ANALYSIS AND DISCUSSION 
3.1. SMBHs 

The primary scientific goal of our Chandra program is to 
measure the luminosities of SMBHs with direct, primary mass 
measurements. For each nuclear source we list the flux and 
luminosity (assuming isotropic radiation and the distance to 
each galaxy listed in Table 1) for the 0.3-2, 2-10, and 0.3-10 
keV bands in Table 3. The fluxes listed are intrinsic (unab- 
sorbed). Note that the fluxes in each band are derived from 
fits to the X-ray spectrum using different bands as normal- 
izations, and, therefore, the uncertainties are correlated. In 
Figure 6 (and 7), we show plots of the spectra of the nuclear 
sources in all 12 galaxies, along with the best-fit model spec- 
trum for each source. The parameters of the best-fit model 
spectra are in Table 2. The spectra have been binned for visu- 
alization purposes. The spectra are all well fit by a power law, 
and none of the spectra requires a more complicated model. 

Figure 2 is a histogram we call a "distributed histogram" of 
the values of the power-law exponent (F) inferred from our 
spectral fits. The uncertainties in F in some spectral fits are 
much larger than some of the others and in some cases much 
larger than the size of the histogram bins in which we are 
interested. To account for this we assume that the errors in 



Table 2 

Results of spectral fits 



ualtixy 


ID 
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M,, 
'V// 


Iogr2_io 
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^"/d.o.t. 




("91 




("41 








NGCI300 


3 




O.Oloo 


-i5.o;g-^ 


2-6;i| 






5 




4.4«-^ 


-i4.5i|!i 


10.0 






6 


ULX 


0.6+°-^ 


-14.6+2-? 

—0.2 


4.3+i-l 

—0.9 






9 


ULX 


o.9+i;| 


-15.5+U-= 


5.9+t; 

—3.4 






10 






-15.0+U-^ 








12 


ULX 


0.3+g-^ 


-15-1+05 


3.8+1-5 






13 


Nuc. 


o-olSi 


-13.0;g;] 


o.4;si 


4.51/3 




14 




—0.6 


-14.8 


4.1 




NGC2748 


18 


ULX 


1-6+0 5 


-i3-2;gl 


2-3+S-? 


2.09/ 3 




20 




o-iioi 


-14-4+0 3 


—0.5 






21 


Nuc. 




-15-0+0 6 


5.0+^-; 






22 


ULX 


l-O+oS 


-i5.o;gi 


5-l+f6 






23 


ULX 


2.3+1-2 


-14.4+g-5 


5-1+" 






24 


ULX 


l-6+o1 


-13-5+gi 


3.0+S-= 






25 


ULX 


0.2«-| 


-13.7;g-| 


2-9;s-j 


6.23/ 9 


NGC2778 


26 


ULX 


0-2+07 


-14.5;g-J 


2-4;i-g 






27 


Nuc. 


0.2+f, 


-15.7+°i 

—0.8 


4.6+?i 

— 1.4 




NGC3384 


29 




Q.Ot^l 


-i3-3;i!i 


i-oiS-? 






30 


ULX 


0-6+03 


-i5-2;gi 


4-1+1-? 






32 


Nuc. 


O-l+oI 


-i3-8;|!i 


2-O+S-J 






34 


ULX 


0.8«-* 


-i4.9;g-= 








35 


ULX 


—0.2 


-15.2+|!-1 

—0.4 


4.2+!-* 

— 1.4 




NGC4291 


36 




o.i+iJ-J 


-14.5«-^ 


1-8+0 6 






37 


ULX 


l-6+oi 

—0.6 


-14.8+2-5 

—0.3 


6.1+?-S 

— 1.7 






38 




0-6+05 


-i4.o;i!i 


1-4+0 9 






39 


Nuc. 


O-l+ol 


-13.9+Si 


1-7+0 J 






40 


ULX 


o.o«-i 


-i3.6;g-| 




1.07/3 
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1 1+5J 

-0.7 


-15 6+2-5 

-0.7 


6 3+" 




NGC4459 


42 


Nuc. 


o-o+S'i 


-13-6+g-l 


1.8+S-? 


4.28/ 4 




43 


ULX 


-0.8 


-14.6+2-? 

—0.2 


4.7+?-; 

-1.7 






44 




0.3+°T 

—0.1 


-13.6+°-| 

—0.1 


1-8+2-t 

—0.3 




NGC4486A 


46 






-i4.2;g-^ 








47 


Nuc. 


3+i;-5 

"--^-0.3 


-14 4+2-i 


^■"-i.o 




NGC4596 


48 


Nuc. 


o-olSs 


-i4.4;g-^ 


2.0+S-2 






49 




o-olSs 


-13-911!-^ 


1-51S1 






50 




0.1+2-! 

-0.1 


-14.0+25 


2.0+»-t 

-0.5 






53 




o.o+°-2 

—0.0 


-13.9+2-; 

—0.2 


1-1+0 3 

-0.3 




NGC4742 


54 




O-SlSi 


-i3.8;g-| 








55 


ULX 


1-1+05 


-14.5+U-^ 


3.5+?-* 






56 


Nuc. 


o.i«i 


-i3.6;g-| 


1-71S1 


7.01/4 




57 


ULX 


0.4«;^ 


-i5-4;gi 


4.3«:5 






58 




n+»-' 




'-^-0.3 




NGC5077 


59 


Nuc. 


n+»-' 




'-^-0.2 


2.18/4 




60 






-15 1+°-^ 


3-1!?:^ 




NGC5576 


61 


ULX 


"-^-0.4 


_14 7+0-5 








62 




4.8 


-14.5 


9.5 






63 


Nuc. 


0-51; 


-14 5+°-^ 


3 5+'-' 






64 


ULX 




-15 7+°-'' 
'-0.8 


5 1+2-* 




NGC7457 


67 


Nuc. 


0-31?:; 


_ 14 4+0.3 


2 5+'-° 






68 




0.3+°-? 


_ 14 5+0.3 

^^-^-0.4 


2.8l„;5 





Note. — This table lists the best-fit parameters with la uncertainties of our 
spectral fits to sources with sufficient counts to waiTant fitting. If we could 
not reliably obtain uncertainties, we omit the listing of errors and consider the 
results approximate. Columns list: (1) the name of the galaxy in which the 
source appears to lie; (2) our running identification number; (3) classification of 
source; (4) total absolution column towards the source in units of 10"^ cm~"; 
(5) logarithmic normalization of the powerlaw in unabsorbed (intrinsic) flux in 
the 2-10 keV band in units of erg s~' cm~"; (6) powerlaw slope; and (7) x~ per 
degrees of freedom if there were sufficient counts to use x~ statistics. 
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Table 3 

Nuclear fluxes 



Gfllflxy 


ID 


f0 3 1 






Fj 10 






ft 3 10 




ioi 10 


^'-'&\^2-10/ '-'had / 


CI 1 


(2) 


(31 




(41 


w> 






(71 




(8) 


(9) 


NGC1300 


13 


7.3;Sj X 10- 


15 


3.S+11 X 10-™ 




13 


5.211? X 10-" 


l.llgixlO- 


13 


5.6+] ? X 10" 


-6.3 ±0.3 


NGC2748 


21 


2.0;'i';t X 10- 


15 


1.9+fg»-^x 10*° 


7.3tf/ X 10- 


16 


6.91JV X 10" 


2.71;;^ X 10- 


15 


2.0+f*-* X 10* 


-8.0±0.7 


NGC2778 


27 


4.3lJ;^ X 10- 


15 


1.9ltVx 10-™ 


1.9lfg%10- 


16 


1.411? X 10" 


4.41J;? X 10- 


15 


1.9l-;VxlO" 


-8.2±0.8 


NGC3384 
NGC4291 


32 
39 


I.3l!J]x 10- 
5.61'j'-^ X 10- 


14 
15 


3.21^* X 10'* 
6.6-^^i X 10'* 


1.7i;jxl0- 
1.4+0-5 X 10- 


14 
14 


2.8l|-^ X 103* 
1.01°^ X lO'' 


3.01;;^ X 10- 
1.912;^ X 10- 


14 
14 


5.91^ * X lO'x 
1.7lgixl03' 


-6.9±0.2 
-7.6±0.4 


NGC4459 
NGC4486A 


42 
47 


i.6;s;jxio- 

I.8;S;*xlO- 


14 
15 


6.3i;-» X 10^* 
1.3lJ'g' X 10'* 


2.3^2^? X 10- 
4.01^;^ X 10- 


14 
15 


8.01^5 X 10'* 
1.4l|!,-^x 10'* 


3.91?;^ X 10- 
5.8l|^ X 10- 


14 
15 


1.4lU-2xl03' 
2.8liV X lO'x 


-7.1 ±0.1 
-7.1 ±0.4 


NGC4596 


48 


4.21?;^ X 10- 


15 


1.8lg-^x 10'* 


3.91^^x10- 


15 


1.5li-Sx 10'* 


8.ll°;J X 10- 


15 


3.31J ? X lO'x 


-7.9±0.3 


NGC4742 


56 




14 


5.2-^?-^ X 10'* 




14 


7.2l[5 X 10'* 


3.8l|;| X 10- 


14 


1.2lg-2xl0» 


-6.4±0.2 


NGC5077 


59 


1.312;^ X 10- 


14 


3.6lg-? X 103' 


2.8^°:? X 10- 


14 


6.8l|-^ X 10-" 


4.11? ^ X 10- 


14 


l.Olg-JxlO* 


-7.2±0.2 


NGC5576 


63 


4.715 5 X 10- 


15 


4.31JV X W 


2.7+51 ^ 


15 


2.6l|-? X 10-« 


7.41J;^ X 10- 


15 


4.61 1;=/ X 1039 


-8.0±0.3 


NGC7457 


67 


3.9i;;{ X 10- 


15 


3.11]-^ X lO'* 


4.3tyf X 10- 


15 


1.0+»|x 10-« 


8.2l^;f X 10- 


15 


4.21^;* X 103* 


-6.7 ±0.3 



Note. — Fluxes and luminosities of nuclear sources, assumed to be the SMBH in each galaxy. For each source, we list the F, the absorbed (apparenti flux 
in units of erg s"^ cm"^, and L, the unabsorbed (intrinsic) luminosity in units of erg s"', for each of the 0.3-2, 2-10, and 0.3-10 bands. Uncertainties are 
listed as 1(7 intervals, and the final column includes uncertainties in the mass. Note that because of covariances between the model parameters, the full band 
is not simply the sum of the soft and hard bands with uncertainties added in quadruture. The final column lists the logarithmic Eddington fraction. While 
some of the sources have full-band luminosities consistent with zero at about the 3(7 level, point sources at the centers of each galaxy at greater than 4(7 
confidence. This is because the luminosity depends upon an unknown spectral form whereas detection depends on raw count rate above the background. For 
example, the full-band luminosity for source 27 in NGC 2778 is only about 1(7 above zero, but inspection of Table 5 shows that the net count rate for this 
source in the full band (summing uncertainties in quadruture) is (9. 6 ± 2.1) x IQ— * ctss"' . 



r are distributed normally and plot the contribution to each 
bin according to the value and error. That is, for each value 
of r, the histogram in a bin of width w that is xa away from 
the central value is w((7\/27r)^' exp(— where a is the 
error of the measurement. We also include analogous data 
from Giiltekin et al. (2009a). There is an obvious peak at 
r « 2. Related to the distributed histogram we introduce a 
"distributed median," which we define as the point at which 
the area under the distributed histogram is half of the total 
area. The distributed median of F is 1.86, very close to the 
canonical 1.7 power-law for AGNs (e.g., Mushotzky 1984). 




-4 -2 2 4 6 



Figure 2. Distributed histogram of values of F, the slope of the power-law. 
This is a histogram of the values of F inferred from our spectral fits, with each 
value spread out among the bins according to the error on the value, assuming 
a normal distribution. That is, for each value of F, the histogram in a bin of 
width n' that is .vfj away from the central value is w 

where C7 is the en'or of the measurement. This method of visualizing the his- 
togram is useful when the errors vary in magnitude from one measurement 
to the next and when the errors are larger than the size of an interesting his- 
togram bin. The new data from this paper are in red, and the new results plus 
those from Giiltekin et al. (2009a) are in gray. The point that divides the area 
under the distributed histogram in two is F = 1 .86, very close to the canonical 
F = 1.7 power-law in AGNs (Mushotzky 1984). 



Because we have direct, primary measurements of the mass 
of the central black hole in each of these galaxies, we can 
report luminosities as true Eddington fractions. In Figure 3, 
we plot a distributed histogram of log(L2-io/i'Edd)- We in- 
clude uncertainties in the black hole mass but instead assume 
the value given in Table 1. As can be seen in the figure, 
the hard X-ray Eddington fractions are small, in the range 
10^** < i2-io/i'Edd < 10^^, and the distributed median of the 
logarithmic Eddington fraction is —7.2. The reason that all 
sources are at such low accretion rates is that the sources 
were selected based on having a dynamical mass measure- 
ment of the central black hole. Most of the black hole masses 
in Giiltekin et al. (2009c) are based on stellar dynamical mod- 
els (e.g., Giiltekin et al. 2009b) and gas dynamical models 
(e.g., Barth et al. 2001), and these methods are best when the 
there is no contamination from AGN light. Typical bolometric 
corrections for low-luminosity sources such as these is ^ 10 
(Vasudevan & Fabian 2007), so that the bolometric Eddington 
fractions of these sources are in the range 10^^-10^^. 

Note that the power of using dynamical mass measurements 
instead of secondary mass estimates is evident in our tabula- 
tion of L2-io/^Edd values. The median error in log(L2-io) is 
0.25 dex, the median error in logMsH is 0. 17 dex, and the me- 
dian error in secondary mass estimates is 0.5 1 dex. Thus the 
uncertainty in secondary mass estimates dominates the total 
uncertainty in log/gdd and is completely subdominant when 
using primary measurements. 

In Figure 4 we plot the spectral power-law slope as a func- 
tion of log(L2-io/i'Edd) for the current data and for the data 
in Giiltekin et al. (2009a). We also fit for a linear relation be- 
tween these two quantities. There is one source, Cen A, with 
L2-io/iEdd > 10^^ and F < 0, that is likely Compton thick 
so that a power-law index at energies below 10 keV possibly 
probes different energetics and/or accretion physics than the 
rest of the sources. Because of this and the potential lever 
arm it could have on the fit, we fit both with and without this 
source, but it had no effect on our results. Since there is a 
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log(L,_,„ / Lem) 

Figure 3. Distributed histogram of values of log(L2-io/^Edd). the slope of 
the power-law. This is a histogram of the values of L2-io/i'Edd inferred from 
our spectral fits, with each value spread out among the bins according to 
the error on the value, assuming a normal distribution. That is, for each 
value of log(L2_io/iEdd) with measurement error a, the histogram in a bin of 
width n' that is xa away from the central value is w{a\/27l)^' exp(— .r^/2). 
This method of visualizing the histogram is useful when the errors vary in 
magnitude from one measurement to the next and when the errors are larger 
than the size of an interesting histogram bin. The distributed histogram shows 
that the accretion rates that we study in this sample are very low. 

Strong covariance between F and L2-10, we do not include 
measurement uncertainties on L2-1Q but instead bootstrap (re- 
sample with replacement) the sample to estimate uncertain- 
ties, expressing our result as the median with 68% interval. 
A histogram of the bootstrap results showed roughly Gaus- 
sian distributions, indicating reliability of results. Finally we 
include a scatter term, assumed normally distributed in F be- 
cause the data clearly deviate from any single line by more 
than their measurement uncertainties. The fitting method and 
code are the same as that used and described in Giiltekin et al. 
(2009c). The method uses a generalized maximum likeli- 
hood method that is capable of including upper limits, arbi- 
trary error distributions, and arbitrary form of intrinsic scatter. 
Here we assume Gaussian errors in F and Gaussian scatter 
in the F direction. We plot the best-fit linear relation, which 

is F = 1.8±0.2- (0.24±0.12)log(L2-io/10"'^LEdd) with an 
rms intrinsic scatter of 0.65 ± 0.20 (uncertainties are 1 (7 here 
and throughout this paper). Our results are inconsistent with 
a slope of zero or larger at the 2(7 level. 

To illustrate what the additional data in this program add as 
well as how sensitive these results are to the use of dynam- 
ical masses, we perform two exercises. First, we repeat the 
fit without the new data, yielding a slope of —0.33 ±0.25. 
So while this is consistent with our full results, it is far less 
conclusive as to whether the X-ray spectral properties at very 
low Eddington rates is different from those at high Eddington 
rates. Second, we fit with the full sample but with masses and 
uncertainties generated from the M-a relation: 

log(MBH/M,,) = a+j3log((7/200kms-i), (3) 

where a is the velocity dispersion of the host galaxy, and 
a = 8. 12 ±0.08 and )3 =4.24 ±0.41 are the best-fit M-a pa- 
rameters (Giiltekin et al. 2009c). Propagating uncertainties. 
Ex, for each quantity, x, the variance in logarithmic mass due 
to random errors is 

eiogMBH = + H (^/200 kms-')]^ej + ^el + ei, (4) 
where Eq = 0.44 is the observed intrinsic scatter in the M-a 



relation. When fitting with these new masses and uncertain- 
ties, the slope is —0.27 ±0.13. So it appears that without 
direct masses, we would have come to the same conclusion, 
though by using masses derived from the M-a relation for 
the very galaxies from which M-a was derived, we necessar- 
ily underestimate the amount of systematic errors introduced 
from using a secondary quantity. 

As mentioned in section 2.1, one consequence of our sam- 
ple selection is that we only cover SMBHs emitting at very 
low Eddington rates. A direct consequence of this is that 
our X-ray observations are more susceptible to contamina- 
tion from a fixed amount of hot gas than would be an X-ray 
source that was intrinsically brighter Since hot gas emission 
typically peaks at an energy below 1 keV, any contamination 
hot gas in our spectral data would tend to make the spectrum 
softer and thus we would infer a larger value for F. This is par- 
ticularly concerning since our results indicate larger values of 
F at lower Eddington rates, as this selection effect and con- 
tamination would manifest. We have mitigated this as much 
as possible with our selection of background regions, which 
are annuli that surround the source regions. This effectively 
removes the contribution of hot gas that is seen in the back- 
ground region. If, however, the diffuse emission is more cen- 
trally concentrated than the background annulus, then there 
will still be some contamination from diffuse gas. To take this 
into account we refit all nuclear sources with an additional as- 
trophysical plasma emission code (APEC Smith et al. 2001) 
component (with abundance fixed to solar). For all but three 
the results for L2-10 and F are consistent at the lc7 or better, 
and the remaining 3 are consistent at about the 2a level. This 
lack of significant change strongly suggests that our results 
are not affected by the contamination of diffuse gas. 
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Figure 4. Slope of the spectral power-law (F) as a function of hard X-ray 
Eddington fraction (L2-lo/iEdd)- Large red circles are new results from tliis 
paper, and small gray squares are results from Giiltekin et al. (2009a). The 
slope of the best-fit relation to these two data sets, which is drawn as a solid 
black line, is 0.24 ±0.12. The dark shaded region shows the lt7 confidence 
band, and the light shaded region shows the 1 <y confidence band plus the rms 
intrinsic scatter Error bars on each point are Iff uncertainties. For reference 
we show similar data based primarily on secondary mass estimates from Gu 
& Cao (2009), Shemmer et al. (2008), Winter et al. (2009), and Younes et al. 
(2011) as indicated in the legend drawn as small open triangles. 

Our results are consistent with previous works looking at 
the anti-correlation between F and Eddington fraction. Us- 
ing a sample of 55 LLAGNs with a mixture of primary and 
secondary black-hole-mass measurements, Gu & Cao (2009) 
fit a linear relation assuming a constant LBoi/i2-io = 30. For 
iBoi/^Edd < 10"^ they inferred F = 1.55 ± 0.07 - (0.09 ± 
0.03)log(LBoi/LEdd)- Note that the differences in intercepts 
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between the fits are all a result of defining the independent 
variable differently. While their results are consistent with 
ours at about the 1.2(7 level, the absence of a scatter term 
in the fitting function may skew the results since the resid- 
uals are larger than would be expected just from measure- 
ment errors. As part of the Chandra Multiwavelength Project, 
Constantin et al. (2009) studied 107 LLAGNs with Chandra 
observations. They found a strong anticorrelation between F 
and LBoi/i-Edd, where they assumed LboI = I6L2-10 and used 
the M-a relation to estimate black hole masses. Their Spear- 
man rank correlation coefficient was —0.75, and their fit to the 
relation was T ^ 0.98 ± 0.13 - (0.27 ± 0.04) log(LBoi/iEdd)- 
The goodness-of-fit, however, was X'^ = 274 for 105 degrees 
of freedom, indicating that it was very unlikely that the data 
could come from a scatter-free model. In a study of 153 
AGNs detected by the Swift Burst Alert Telescope hard X- 
ray instrument. Winter et al. (2009) found a no correlation 
between F and L2-io/i'Edd- It is possible that their sample, 
which is primarily at L2-io/i'Edd > lO^'* and reaches up to 
L2-io/i'Edd ~ 0.04 is actually measuring an energetically dif- 
ferent mode of accretion than our sample. In any case, our 
study agrees with Winter et al. (2009) that the spectra do not 
harden with increasing accretion rates. In a study of 13 LIN- 
ERs with XMM-Newton and/or Chandra observations and as- 
suming very different M-a relation (due to Graham et al. 
201 1), Younes et al. (201 1) found F = 0.1 1 ±0.40 - (0.31 ± 
0.06)log(L2-io/i'Edd)- A similar anti -correlation is suggested 
in X-ray binary data (Corbel et al. 2006). Given the differ- 
ence in assumptions and fitting techniques, we consider all 
of these results consistent with each other in the following 
conclusions: (1) there is an anti-correlation between the hard 
X-ray photon index, F, and Eddington fraction for LLAGNs, 
and (2) the correlation between F and log(L2-io/i'Edd) has a 
slope of roughly —0.2. 

The anticorrelation between F and L2-io/iEdd is in contrast 
to AGNs emitting at higher Eddington fractions. For example, 
Shemmer et al. (2008) found from fits to a sample of 35 radio 
quiet, luminous, and high-Eddington-fraction (Leoi/i-Edd ^ 
0.01-1) AGNs that there was a strong positive correlation. 
From their fits, they found log(LBoi/i'Edd) = —2.4 ±0.6 + 
(0.9 ± 0.3)F, corresponding to F = 2.6 + l.l(LBoi/i'Edd)- 
Thus our evidence leads us to conclude that we are seeing 
a change in the physical processes responsible for emission at 
low Eddington rates. 

The softening of the spectrum with decreasing Eddington 
fractions at low mass accretion rates as we find in our sample 
is predicted by advection dominated accretion flow (ADAF) 
models in stellar-mass X-ray binaries (Esin et al. 1997), which 
should reasonably translate to the low accretion rates seen 
in our SMBH sources. Esin et al. (1997) models predict a 
steepening of the 1-10 keV band index of F « 1.7 to about 
2.2 over mass accretion rates of log(m) = —2.5 to —4, or 
a slope of approximately 0.33. The comparison is not di- 
rect, but our results are fully consistent with this prediction, 
which is attributed to the fact that at lower accretion rates 
bremsstrahlung emission becomes relatively more important 
compared to Comptonizaion. 

3.2. ULXs 

Ultraluminous X-ray sources (ULXs) are a category of X- 
ray emitting point sources that are too bright to be explained 
by isotropic emission resulting from sub-Eddington accretion 
onto stellar-mass (M < 10 Mq) black holes and are also non- 



Table 4 

Ultraluminous X-ray Source Candidates 



Galaxy 


ID 




P(L>LvLx) 


if 


p 


(Wbg) 


(1) 


(2) 


(3) 


(4) 


(5) 


(5) 


(6) 


NGC1300 


6 


8.91^0° X 10" 


0.92 


X 10-'* 


0.07 


0.13 




9 


2.8;^''g"i ' X 10"" 


0.50 


2.31S-? X 103* 


0.83 


0.13 




12 




0.43 


3.51°-* X 103* 


0.74 


0.13 


NGC2748 


18 


1.4+^ ? X 10"" 


>0.99 


7.0+i-;, X 103' 


0.002 


0.02 




22 


2.8+™-' X 10™ 


0.88 


5.2«-S X 103* 


0.05 


0.02 




23 


9.2;J7" X 10™ 


0.91 


73tll X 103* 


0.03 


0.02 




24 


X 10" 


>0.99 


4.6;S * X 103' 


0.002 


0.02 




25 


l.Otll X 10™ 


>0.99 


5.S+_°at X 103' 


0.002 


0.02 


NGC2778 


26 


8.0+^V X W 


0.16 


63+_ji X 103* 


0.79 


0.01 


NGC3384 


30 


51+";^ X lO's 


0.31 


e.otll X 103' 


0.16 


0.03 




34 


5.6+™*" X lO'* 


0.37 


9.4+5 J X 103' 


0.39 


0.03 




35 


7.0+f/ X W 


0.33 


1.3;»-3x 103* 


0.49 


0.03 


NGC4291 


37 


2.91^7' X 10"' 


0.99 


7.8;|-^ X 103* 


0.008 


0.01 




40 


l.&tTi X lO'" 


0.96 


2.5;S-^ X 103' 


0.60 


0.01 


NGC4459 


43 


1.7™ X 10™ 


0.77 


3.3;[-3x 103* 


0.03 


0.01 


NGC4596 


52 


1.5if5' X 10"' 


0.85 


2.4+2'^ X 103* 


0.27 


0.01 


NGC4742 


55 


1.6+fj»x 103" 


0.46 


2.5+i-^ X 103* 


0.16 


0.00 




57 


9.6;«^-3 X 103« 


0.41 


1.3lS|x 103* 


0.72 


0.00 


NGC5576 


61 


1.3+^'g''x 103" 


0.42 


5.7+^« X 103* 


0.57 


0.03 




64 


5.9;™ ' X \0"> 


0.63 


4.5t?-? X 103* 


0.70 


0.03 



Note. — A listing of the ULX candidates identified in this survey. Columns list: (1) the 
name of the galaxy in which the source appears to lie; (2) our running identification number; 
(3) the unabsorbed (intrinsic) luminosity in the 0.3-10 keV band; (4) the probability, based 
on the flux parameter uncertainty, that the given source's luminosity is above the definition 
of a ULX (Lo,3-lo > 2 x 103' erg s"'); (5) the 0.3-10 keV luminosity inferred when fixing 
the absoiption column to the Galactic value towards that source; (6) the p-value result of 
our simulations to calculate the significance of the improvement of including as a free 
parameter; and (7) the expected number of background sources in the galaxy with fluxes 
^0,5-2 >2x 103' erg s"' (4;rZ)")"' . The probability of having at least one background source 
in the galaxy of such flux is 1 — exp( — (Wbg)) ^ (A'bg)- All uncertainties are listed as IfT 
intervals. 

nuclear so that they are unlikely to be the galaxy's central 
SMBH. ULXs are an intriguing class of source because if 
(1) they are emitting roughly isotropically (i.e., not strongly 
beamed towards our line of sight), so that we may correctly in- 
fer their luminosity, and (2) they are not accreting well above 
the Eddington limit, so that we may robustly infer a lower 
limit to the mass, then the most natural explanation is a black 
hole of mass ^ 10^-10^ M. . These intermediate-mass black 
holes (IMBHs) would fill the gap in mass between stellar- 
mass black holes and SMBHs. IMBHs are interesting because 
their formation in the local universe requires a non-standard 
path (Miller & Hamilton 2002; Giiltekin et al. 2004, 2006; 
Portegies Zwart et al. 2004). There are alternative interpre- 
tations to ULXs other than IMBHs, including beamed, non- 
isotropic emission (King et al. 2001). This would change the 
luminosity inferred from the flux so that it is consistent with 
sub-Eddington accretion onto stellar-mass black holes. There 
are at least a few sources (Kaaret et al. 2004; Pakull & Mirioni 
2001) where emission from the surrounding medium argues 
in favor of roughly isotropic emission. Another alternative is 
to have super-Eddington accretion, which has been invoked 
in a number of different ways (Begelman 2002, 2006). The 
very bright source ESO 243-49 HLX-1 is similarly difficult to 
interpret without invoking an IMBH (Farrell et al. 2009). 

Although our survey was not targeted at ULXs, we are sen- 
sitive to them and present a list of ULX candidates in table 
4. We adopt the definition of (Irwin et al. 2003) that ULXs 
are sources with L0.3-10 > Lulx = 2 x 10^^ erg s^', which 
avoids contamination from bright, massive stellar-mass X-ray 
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binaries. We assume that all sources are isotropically emit- 
ting at the distance of the host galaxy. Because of measure- 
ment uncertainties, there is always a finite probability that a 
source that appears to have Lo.3-10 > ^ulx, is actually in- 
trinsically too dim to be a ULX. Additionally, sources that 
are just below Lulx still viable ULX candidates. For 
these reasons, we list all sources that are at least 1 a consis- 
tent with being a ULX and list P(Lo.3-io > ^^ulx), the prob- 
ability of the source having L0.3-10 > i'ULX- This is calcu- 
lated by finding AC, the change in fit statistic, when setting 
■Pb.3-10 = 2 X lO"*^ erg s"' (47rD^)"' and refitting. Then we 



calculate P(Lo.3-io > ^ulx) = 0.5 ± 0.5 erf (AC/2) 



nO.5 



, tak- 



ing the top and bottom sign for when the best-fit luminosity is 
above or below Lulx, respectively. 

One source of contamination is background AGNs that ap- 
pear to be in the galaxy. We have checked the catalog of 
Veron-Cetty & Veron (2010) and found no known AGNs at 
the locations of our ULX candidates, but it is clearly possible 
for a previously unknown AGN to be located at this position. 
To quantify this effect, we calculate {Nqg), the expected num- 
ber of background AGNs in the region of the galaxy target 
with fluxes greater than 2 X 10^^ erg s^' (4;rD^)^'. Todothis 
we use the fits to the background AGN density as a function 
of flux from (Giacconi et al. 2001). Unfortunately, they only 
provide fits for fluxes in the 0.5-2 and 2-10 keV bands. Since 
most of the flux in a power-law source comes from the soft 
band, we use the softer band and assume that all of the flux 
comes from this band. Since the relations are approximately 
linear, if only half of the emission comes from this portion 
of the band, then the average number of background sources 
would roughly double, still a small number. It is important 
to use the whole area of the galaxy and the smallest possi- 
ble flux for the calculation of (A^bg) because we would have 
listed any source above the threshold flux apparently within 
the galaxy as a ULX candidate. The expected background for 
each galaxy is listed in Table 4, and it is small. The proba- 
bility of a galaxy having at least one confusing background 
source is 1 — exp(— (A^bg)), and a lower limit f^or the proba- 
bility of an individual ULX candidate's truly being a ULX is 

P(ULX) -P(Lo.3-io >iuLx)exp(-(A?BG)). 

Note that the probability of ULX luminosity that we calcu- 
late is only accurate if our spectral model is a reasonably good 
model. In Figure 8 we plot spectra of the six sources (6, 18, 
23, 24, 25, and 40) with P(Lo.3-io > iuLx) > 0.9. Sources 18, 
23, and 24 all have A^h > 10^^ cm^^ and thus while their ap- 
parent flux is rather modest, the inferred, intrinsic absorption- 
corrected luminosity is quite high. This inference depends 
strongly on the correct estimation of A^h, which derives from 
the assumption of a power-law spectral form. We have tried 
fitting with other spectral forms for these three sources but 
were not able to come up with acceptable fits. In particular, 
we tried more complicated models with combined disk black- 
body and power-law spectral forms, but the disk blackbody 
temperature normalization was either unphysically high or its 
normalization was so low as to make the component irrele- 
vant. 

To address the issue of a ULX classification's sensitivity to 
the unknown intrinsic absorption column, we refit all ULX 
candidate spectra with A^h fixed to the Galactic value in Table 
1 for each galaxy. The resulting luminosities (Lp) are listed 
in Table 4. We also calculated /^-values of a likelihood ratio 
statistic using Monte Carlo simulations using the method of 



Protassov et al. (2002). Each simulation synthesized 1000 re- 
alizations of the ULX spectra based on the best-fit fixed-A^H 
model. We then fitted each synthetic data set with a fixed- 
A^H model and one in which A^h was a free parameter, and 
calculated a likelihood ratio for each spectrum. This proce- 
dure generated a distribution of our likelihood ratio statistic 
according to the null hypothesis that A^h fixed at the Galactic 
value is the correct model. We adopt p < 0.01 as our level of 
significance for requiring an additional free parameter Using 
the fixed A^h spectral model, only 4 sources (18, 24, 25, and 
40) are bright enough to be considered ULXs. We note that 
source 25 has a luminosity of L0.3-10 = 1.03 x lO'*'^ erg s^', 
with P(Lo.3-io > iuLx) > 0.99 and is very well fit by the ab- 
sorbed power-law model (p < 0.002). 

Finally, we note that NGC 4291 is also represented in the 
XMM-Newton catalog of ULXs due to Walton et al. (2011). 
In NGC 4291 we find 2 ULX candidates (sources 37 and 39). 
Given the best-fit absorption column, both of these are likely 
to have a flux bright enough to be true ULXs (1(Lo.3-io > 
iuLx) > 0.96), and source 40 has Lp = 2.5^11 x 10^*^. Nei- 
ther source, however, is listed as a ULX in Walton et al. 
(201 1). These two sources are within 9"of the nuclear source 
and would not be reliably resolved with XMM-Newton, and 
Walton et al. (2011) did not consider bright sources within 
15" of the center of elliptical galaxies or any sources within 
7"5 of the center of the galaxy. On the other hand, Walton 
et al. (2011) fist 2XMM J122012.5+752204 as a ULX can- 
didate, using a definition of Lulx = 1 x 10^^ erg s^^ This 
XMM source is consistent with the position of our source 36 
(CXOU J 122012. 1+752203) for which we measure a lumi- 
nosity Lo.3-10 = 4.0+j g X 10^** erg s"'. The Walton et al. 
(2011) classification of this source as a ULX is based on the 
2XMMS serendipitous source catalog (Watson et al. 2009) 
0.2-12 keV flux of (5.5 ± 1.0) x lO""* erg s^^ cnT^, which, 
at our adopted distance, corresponds to an isotropic lumi- 
nosity of (4.1 ±0.8) X lO^**. The 2XMMS flux measure- 
ment assumes a spectral form of an absorbed power-law with 
r = 1.7 and A^h = 3 x 10^*^ cm"^, very close to the best-fit 
values we find from our spectral fit. If we fix F = 1 .7 and 
A^H = 3 X 10^" cm^^ and fit to our Chandra data, we get a 
0.2-12 keV flux of (4.9± 1.5) x lO^'^ erg s^^ cm^^. Given 
that ten years elapsed between the XMM observation (MJD 
51671) and the Chandra observation (MJD 55541), we con- 
clude that the source is variable on these scales. 

3.3. Off-nuclear sources 

In addition to the nuclear and ULX candidate sources, we 
detected 36 off-nuclear sources whose position on the sky is 
consistent with being in the galaxy. As the flux goes down, the 
probability for a source to be a background AGN increases. 
Thus we cannot be certain that the sources are intrinsic to the 
galaxy, but, generally, the clustering of point sources within 
the galaxies' optical extent is much higher than over the en- 
tire Chandra image. To summarize the off-nuclear sources, 
in Figure 5 we plot Co.3_i/Ci_2, the ratio of the count rates 
in the 0.3-1 to 1-2 keV bands, as a function of C2-10, the 2- 
10 keV count rate. With only 2 data points (one color and 
one intensity), it is not possible to break the degeneracy in an 
absorbed power-law model between A^h and F, but the range 
of parameters needed to reproduce most of the data points is 
reasonable. 

4. SUMMARY 
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Figure 5. Color-intensity plot of off-nuclear sources. As a function of C2-10, 
the 2-10 keV count rate, we plot the ratio of the count rates in the 0.3-1 to 
1-2 keV bands. Limits are given with arrows at the 3f7 limit. We do not 
plot source 4, which only had upper limits in all three bands. The horizontal 
lines are color-intensity curves for an absorbed power-law model with the 
parameters and indicated in the legend (A^h in units of 10"^ cm^^). With 
only 2 data points (one color and one intensity), it is not possible to break 
the degeneracy between Wh and F, but the range of parameters needed to 
reproduce most of the data points is reasonable. 

In this paper we have presented an X-ray survey of 12 
galaxies with central black hole mass measurements. The ob- 
servations were designed to characterize the nuclear source of 
each galaxy but were sensitive to much more. 

1. Each galaxy was observed for 30 ksec with Chandra, 
and in total we detected 68 point sources in the region 
of the sky occupied by these galaxies. 

2. We detected all 12 nuclear sources with sufficient count 
rates to model their spectra and determine their X-ray 
luminosities and Eddington ratios. The sources all were 
found to be emitting in the 2-10 keV band at 10^**- 
10"*' of Eddington. 

3. When fitting with an absorbed power-law spectral 
model, we found F, the photon spectral index. Fitting 
for F as a function of Eddington fraction, we found 
a negative correlation, consistent with several earlier 
reports on low luminosity AGNs. Our best fit was 
F = 1 .8 ± 0.2 - (0.24± 0. 12) log(L2-io / lO^^LEdd) with 
an rms intrinsic scatter of 0.65 ± 0.20, which is consis- 
tent with predictions from ADAF models, which expect 
bremsstrahlung emission to become more important at 
lower accretion rates. 

4. Our observations were also sensitive to ULXs in the 
target galaxies. We found 20 ULX candidates. Based 
on considerations of the probability distribution of their 
intrinsic fluxes and the probability of having a back- 
ground AGN of sufficient brightness to appear as a 
ULX, we concluded that 6 of these candidates are likely 
(> 90% chance) to be true ULXs. The most promising 
ULX candidate is in NGC 2748 and has an isotropic 
luminosity of L0.3-10 = 1 -03^0:27 x 10"^° erg s"'. 

5. We also present a color-intensity plot of the remaining 
point sources, most of which are likely to be X-ray bi- 
naries local to the galaxy. 

This work will be followed up with EVLA observations of 
the nuclear sources. When combining radio and X-ray data. 



we will be able to provide a complete mass-calibrated funda- 
mental plane that will allow for the estimation of black hole 
masses using X-ray and radio observations. We will also, with 
the possible addition of archival sub-millimeter data, broad- 
band spectral energy distribution modeling. 
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Figure 6. X-ray spectra, folded through the instrument response, of the first half of the sample. The spectra are binned for visualization purposes her The solid 
line shows the best-fit inodel. The second half of the sample is in Figure 7. Error bars show Iff uncertainties. 
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Figure 7. Same as Figure 6 for the second half of the sample. Error bars show la uncertainties. 
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Figure 8. X-ray spectra of strong ULX candidates. Error bars show lc7 uncertainties. The solid lines are power-law models in whicli the Nh parameter was 
allowed to vary above the Galactic value; the dashed lines are models with Nh fixed at the Galactic value towards each source. 
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Table 5 

Chandra point source detections 



Galaxy 


ID 


Source name 


RA 


Dec 


Class. 


0.3-1 keV 


1-2 keV 


2-10 keV 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGCI300 


1 


CXOU J03 1935.0- 192509 


03:19:35.07 


-19:25:09.9 




4.5+1.3 


0.7+0.6 


<3.2 




2 


CXOU J031935. 1-192417 


03:19:35.10 


-19:24:17.4 




0.8+0.6 


1.5+0.8 


<3.2 




3 


CXOU J031935.4-192417 


03:19:35.49 


-19:24:17.6 




2.0+0.8 


2.7+1.0 


0.5+0.5 




4 


CXOU J03 1 936.7 - 1 92500 


03:19:36.78 


-19:25:00.3 




<2.0 


<2.0 


<3.8 




5 


CXOU J031937.8-192607 


03:19:37.88 


-19:26:07.8 




0.9+0.6 


5.6+1.4 


6.0+1.6 




6 


CXOU J031937.9-192441 


03:19:37.95 


-19:24:41.4 


ULX 


4.7+1.3 


11.0+1.9 


2.4+1.1 




7 


CXOU J031938.2-192413 


03:19:38.21 


-19:24:13.5 




1.3+0.7 


2.0+0.8 


0.5+0.5 




8 


CXOU J031938.4-192318 


03:19:38.45 


-19:23:18.0 




1.6+0.8 


2.4+1.0 


<3.7 




9 


CXOUJ031938.5-192452 


03:19:38.58 


-19:24:52.4 


ULX 


2.5±1.0 


4.0±1.2 


1.8±1.0 




10 


CXOU J031939.3-192422 


03:19:39.37 


-19:24:22.8 




2.2±0.9 


3.2±1.1 


2.9±1.1 




11 


CXOU J031939.4-192616 


03:19:39.47 


-19:26:16.0 




<2.6 


1.6±0.8 


3.2±1.2 




12 


CXOU J03 1939.9- 192550 


03:19:39.93 


-19:25:50.4 


ULX 


3.2±1.1 


5.8±1.4 


2.9±1.2 




13 


CXOU J031941.0-192440 


03:19:41.06 


-19:24:40.3 


Nuc. 


5.5±1.5 


12.0±2.1 


27.3±3.1 




14 


CXOU J031941.3-192334 


03:19:41.37 


-19:23:34.2 




3.0±1.0 


5.7±1.4 


1.2±0.8 




15 


CXOU J031942.5-192450 


03:19:42.50 


-19:24:50.7 




0.9±0.7 


3.7±1.2 


3.0±1.1 




16 


CXOU J031942.9-192511 


03:19:42.95 


-19:25:11.2 




1.2±0.7 


1.7±0.8 


3.3+1.3 




17 


CXOU J03 1 945 .9 - 1 92449 


03:19:45.98 


-19:24:49.3 




<2.4 


3.1 + 1.1 


<4.1 


NGC2748 


18 


CXOU J09 1 337.4+7628 1 1 


09:13:37.45 


+76:28:11.4 


ULX 


0.6+0.5 


19.8+2.6 


24.3+2.9 




19 


CXOU J091339.6+762800 


09:13:39.64 


+76:28:00.5 




< 2.6 


1.9+0.8 


3.0+1.1 




20 


CXOU J09 1 340.0+7628 1 3 


09:13:40.01 


+76:28:13.3 




1.9+0.8 


3.6+1.1 


2.0+0.9 




21 


CXOU J091343.0+762831 


09:13:43.04 


+76:28:31.8 


Nuc. 


1.9+0.9 


4.5+1.3 


1.1+0.8 




22 


CXOU J091343. 8+762835 


09:13:43.80 


+76:28:35.9 


ULX 


1.5+0.8 


4.9+1.3 


1.7+0.8 




23 


CXOU J09 1344.4+762829 


09:13:44.42 


+76:28:29.3 


ULX 


0.6+0.5 


4.6+1.3 


3.4+1.1 




24 


CXOU J091348.5+762902 


09:13:48.52 


+76:29:02.2 


ULX 


0.3±0.4 


15.1±2.3 


15.0±2.3 




25 


CXOU J091348.9+762828 


09:13:48.99 


+76:28:28.5 


ULX 


30.3±3.2 


44.4±3.9 


16.7±2.4 


NGC2778 


26 


CXOU J091222.4+350135 


09:12:22.45 


+35:01:35.2 


ULX 


3.7±1.1 


4.7±1.3 


3.3±1.1 




27 


CXOU J091224.4+350139 


09:12:24.40 


+35:01:39.4 


Nuc. 


4.8±1.4 


3.1±1.2 


1.7±1.0 


NGC3384 


28 


CXOU J104813.9+123839 


10:48:13.94 


+12:38:39.4 




0.6±0.5 


3.4±1.1 


2.3±1.0 




29 


CXOU J104815.1+123658 


10:48:15.18 


+12:36:58.7 




6.1±1.5 


16.6±2.4 


12.7±2.2 




30 


CXOU J104816.1+123658 


10:48:16.15 


+12:36:58.7 


ULX 


1.4±0.7 


2.8±1.0 


<2.6 




31 


CXOU J104816.2+123644 


10:48:16.27 


+12:36:44.1 




0.7±0.5 


2.1+0.9 


<2.0 




32 


CXOU J 1048 16.9+ 123745 


10:48:16.97 


+ 12:37:45.7 


Nuc. 


10.7+2.0 


18.7+2.6 


9.8+2.0 




33 


CXOU J 1 048 1 7.4+ 1 23720 


10:48:17.47 


+ 12:37:20.8 




1.0+0.6 


2.0+0.9 


1 .0+0.7 




34 


CXOU J 1 048 1 9.4+ 1 23646 


10:48:19.41 


+ 12:36:46.0 


ULX 


1.0+0.6 


2.4+0.9 


3.6+1.2 




35 


CXOU J104820.8+123839 


10:48:20.81 


+ 12:38:39.4 


ULX 


3.8+1.2 


4.5+1.3 


0.7+0.6 


NGC4291 


36 


CXOU J122012.1+752203 


12:20:12.16 


+75:22:03.1 




1.3+0.8 


1.5+0.8 


1.8+0.8 




37 


CXOU J122016.1+752218 


12:20:16.12 


+75:22:18.0 


ULX 


2.1 + 1.3 


7.7+1.8 


1.1+0.9 




38 


CXOU J122016.5+752151 


12:20:16.50 


+75:21:51.3 




< 3.6 


2.1+0.9 


4.3+1.3 




39 


CXOU J122017.8+752214 


12:20:17.84 


+75:22:14.8 


Nuc. 


6.6+1.7 


7.7+1.7 


6.6+1.5 




40 


CXOU J122019.8+752216 


12:20:19.83 


+75:22:16.2 


ULX 


9.3±2.4 


15.0±2.5 


10.4±2.0 




41 


CXOU J122023.6+752202 


12:20:23.67 


+75:22:02.5 




1.0±0.9 


3.4±1.2 


<3.2 


NGC4459 


42 


CXOU J122900.0+135842 


12:29:00.02 


+13:58:42.0 


Nuc. 


18.5±2.6 


17.2±2.5 


13.6±2.3 




43 


CXOU J122900.5+135825 


12:29:00.56 


+13:58:25.9 


ULX 


1.2±0.7 


4.3±1.2 


2.5±1.0 




44 


CXOU J122900.6+135838 


12:29:00.62 


+13:58:38.3 




4.8±1.3 


12.6±2.1 


10.9±1.9 




45 


CXOU J122900.8+135843 


12:29:00.88 


+13:58:43.6 




<3.4 


2.4±1.0 


<2.4 


NGC4486A 


46 


CXOU J123057.7+121616 


12:30:57.77 


+12:16:16.3 




3.5±1.3 


5.5±1.6 


3.3±1.2 




47 


CXOU J123057.8+121614 


12:30:57.87 


+12:16:14.5 


Nuc. 


2.4±1.0 


2.6±1.1 


3.2±1.1 


NGC4596 


48 


CXOUJ123955.9+101033 


12:39:55.99 


+10:10:33.7 


Nuc. 


4.4±1.3 


4.1+1.2 


2.2+0.9 




49 


CXOU J123956.0+101036 


12:39:56.05 


+ 10:10:36.0 




4.0+1.2 


7.3+1.5 


4.1 + 1.2 




50 


CXOUJ123956.4+101026 


12:39:56.43 


+ 10:10:26.0 




6.0+1.5 


9.5+1.8 


5.1 + 1.3 




51 


CXOU J123956.7+101101 


12:39:56.72 


+ 10:11:01.6 




0.5+0.5 


1.2+0.7 


1.9+0.9 




52 


CXOU J123956.8+101053 


12:39:56.83 


+ 10:10:53.1 


ULX 


4.7+1.3 


2.6+1.0 


1.8+0.9 




53 


CXOU J123956.9+101029 


12:39:56.90 


+ 10:10:29.8 




3.4+1.1 


3.5+1.1 


5.7+1.4 


NGC4742 


54 


CXOU J125147.6-102722 


12:51:47.60 


-10:27:22.6 




4.9+1.3 


12.1 + 1.9 


8.7+1.7 




55 


CXOU J125147.9-102719 


12:51:47.92 


-10:27:19.7 


ULX 


0.6+0.5 


4.3+1.2 


3.2+1.1 




56 


CXOU J125148.0- 102717 


12:51:48.07 


-10:27:17.2 


Nuc. 


14.9±2.2 


18.0±2.4 


12.9±2.1 




57 


CXOU J125148.2-102710 


12:51:48.27 


-10:27:10.8 


ULX 


2.4±0.9 


3.6±1.1 


1.6±0.9 




58 


CXOU J125149.3-102727 


12:51:49.34 


-10:27:27.4 




8.4±1.6 


15.5±2.2 


13.3±2.1 


NGC5077 


59 


CXOUJ131931.6-123925 


13:19:31.66 


-12:39:25.1 


Nuc. 


12.1±2.5 


16.4±2.6 


13.6±2.3 




60 


CXOUJ131931.9-123938 


13:19:31.93 


-12:39:38.3 




1.7±0.9 


2.1±0.9 


1.6+0.9 


NGC5576 


61 


CXOU J142102.7+031611 


14:21:02.76 


+03:16:11.6 


ULX 


0.8±0.6 


3.2±1.1 


2.9±1.1 




62 


CXOU J142102.9+031621 


14:21:02.97 


+03:16:21.6 




1.0±0.7 


3.5±1.2 


1.1+0.8 




63 


CXOU J142103.7+031614 


14:21:03.71 


+03:16:14.9 


Nuc. 


5.1±1.4 


8.4±1.8 


2.7±1.1 




64 


CXOU J142105.1+031615 


14:21:05.18 


+03:16:15.9 


ULX 


2.6±1.0 


2.6±1.0 


2.5±1.0 
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Table 5 — Continued 



Galaxy 


ID 


Source name 


RA 


Dec 


Class. 


0.3-1 keV 


1-2 keV 


2-10 keV 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC7457 


65 


CXOU J230058.0+300850 


23:00:58.07 


+30:08:50.5 




0.7+0.5 


3.1 + 1.1 


2.8+1.1 




66 


CXOU J23()059. 1 +300907 


23:00:59.16 


+30:09:07.4 




0.6+0.5 


1.7+0.8 


6.5+1.6 




67 


CXOU J230059.9+300841 


23:00:59.95 


+30:08:41.8 


Nuc. 


3.2+1.1 


7.5+1.6 


4.3+1.4 




68 


CXOU J230101. 1+300900 


23:01:01.19 


+30:09:00.7 




3.4+1.1 


5.2+1.4 


3.6+1.2 



Note. — List of X-ray point sources detected in the field of each galaxy. Columns list: (1) the name of the galaxy in which the sources appear to lie; (2) a running identification 
number used in this paper; (3) lAU approved source name for each source; (4) and (5) J2()00 coordinates for each source; (6) classification of each source, where Nuc. indicates the 
source is nuclear source and presumed to be the SMBH for its galaxy, and ULX indicates that it is a ULX candidate; (7), (8), and (9) are the count rates for each source in the given 
bands in units of 10""* cts s~^. If a source is consistent with no flux in a given band, then we list the 3<7 upper limit, otherwise we list the 1(7 uncertainty. 



